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We extract a numerical value for the strong coupling con- 
stant a s from the r-lepton decay rate into nonstrange par- 
ticles. A new feature of our procedure is the explicit use of 
renormalization scheme invariance in analytical form in order 
to perform the actual analysis in a particular renormalization 
scheme. For the reference coupling constant in the MS-scheme 
we obtain a s (M T ) — 0.3184 ±0.0060 ea :p which corresponds to 
as(Mz) = 0.1184 ± 0.0007^ ± 0.0006^ mass- This new nu- 
merical value is smaller than the standard value from r-data 
quoted in the literature and is closer to a 3 (A/z)-values ob- 
tained from high energy experiments. 

ll.10.Hi, 12.38.-t, 13.35. Dx 



The physics of r-lcpton hadronic decays is an impor- 
tant area of particle phenomenology where the theory of 
strong interaction (QCD) can be confronted with exper- 
iment to a very high precision. The central quantity of 
interest in this process is the spectral density of hadronic 
states related to the two-point correlator of hadronic cur- 
rents with well established and simple analytic proper- 
ties. The accuracy of experimental data for a variety of 
observables of the r-lepton system is rather good and 
is steadily improving The spectral density itself 

(more precisely, the two-point correlator of hadronic cur- 
rents in the Euclidean domain) has been calculated with 
a very high degree of accuracy within perturbation theory 
(see e.g. Q-Q]). Nonperturbative corrections to the corre- 
lator are known to be small and under control within the 
operator product expansion and factorization approxima- 
tion |]||. The observables in the r system are inclusive 
in nature which makes the comparison of experimental 
data with theoretical calculations very clean |i7|-|i"6|] . Of 
some particular interest is the precise determination of 
the numerical value of the strong coupling constant at 
the low energy scale of the r-lepton mass. Within the 
renormalization group approach this number can then be 
evolved to high energies. This is a powerful consistency 
check of QCD since one is comparing hadron physics at 
a tremendous variety of scales, from one to hundreds of 
GeV (e.g. @). 

In the present note we provide a thorough analysis of 
the procedure of extracting numerical values of a s from 
r-data in perturbation theory. On the theory side one 
expects a high degree of accuracy in the determination 
of a s because of the existence of very accurate pertur- 



bation theory formulas and the simplicity of the renor- 
malization group treatment of the massless quark case. 
However, the numerical value of the expansion parame- 
ter a s is not small at the M T scale and the contribution 
of higher order terms in the perturbation theory series 
can be significant. Arguments have been brought forth 
that the accuracy of finite-order perturbation theory is 
already close to its asymptotic limit which makes the 
interpretation (usually called resummation) of the per- 
turbation theory series in higher orders necessary [p.8| . 
The resummation of contributions related to the running 
of the coupling constant is most advanced e.g. Jl9|-|25f|. 

The decisive new point of our analysis is the explicit 
use of renormalization group invariance in the analysis 
of the r-lepton decay rate within perturbation theory. 
Renormalization group invariance is a fundamental prop- 
erty of perturbation theory in quantum field theory which 
is related to the freedom in defining the subtraction pro- 
cedure f26| . It should be respected in any numerical 
analysis. Renormalization group invariance allows one 
to formally perform the numerical analysis in any renor- 
malization scheme because all schemes are connected by 
a renormalization group transformation. However, in the 
finite-order perturbation theory approach this equiva- 
lence is only approximate due to the systematic omis- 
sion of higher order terms in the perturbation theory 
expressions. This inroduces numerical differences into 
the results obtained in different renormalization schemes. 
Generally one can consider two ways of using perturba- 
tion theory calculations. One is to find relations between 
physical observables which are renormalization group in- 
variant. Then perturbation theory calculations are just 
a purely intermediate step for finding relations between 
observables (see, e.g. j^,^7|) and no numerical analy- 
sis for renormalization scheme noninvariant quantities is 
performed. Indeed, let the perturbation theory expres- 
sions for two observables Oi^ in a given scheme have the 
form 



01 = a s +r ia 2 s +0(a 3 s ), 

2 = a s + r 2 a 2 s + 0(a 3 e ) . 



(1) 



Then the perturbation theory relation between observ- 
ables 0i 2 reads 



2 = 0i + (r 2 - n)0 2 + 0(0\) 



(2) 



and is scheme-independent. The difference r 2 — T\ takes 
the same value for calculations in any scheme. Another 
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way of using perturbation theory calculations is to ex- 
tract numerical values for renormalization scheme non- 
invariant quantities (as the coupling constant in a fixed 
scheme). These are then compared with the results of 
other experiments. In this case the truncation of the 
perturbation theory series leads to numerical violations 
of renormalization scheme invariance and plays an es- 
sential role. In our simple example this means that the 
relations in eq. ([I]) are treated as quadratic functions of 
a a in some fixed scheme and the accuracy of extraction 
of the coupling constant value (and prediction of other 
observables) depends drastically on the scheme used, i.e. 
on the numerical values of the coefficients r\^. 

In the present paper we consider just this second ap- 
plication and extract a numerical value for the coupling 
constant which is not an immediate physical quantity. By 
convention the reference value of the coupling constant 
that is used to compare between different experiments 
is fixed to be the MS-scheme one. However, and this is 
our point in this paper, this does not necessarily mean 
that for its extraction from a given experiment the nu- 
merical analysis should be performed in the MS-scheme. 
It can be more convenient (and numerically accurate) to 
analyze the system in its internal scheme and after find- 
ing numerical values for the internal parameters translate 
them into the MS-scheme using renormalization scheme 
transformation. This program heavily uses explicit renor- 
malization scheme covariance of the theory. However, ex- 
pressions for the amplitudes are available only in pertur- 
bation theory as a truncated series in the coupling con- 
stant. For a truncated series the renormalization scheme 
invariance is only approximate with a precision of the 
order of the value of the first omitted term. Therefore 
numerical values obtained in the MS-scheme directly and 
through renormalization group transformations can dif- 
fer. We discuss this problem and argue that the inter- 
nal scheme results are most reliable physically and are 
more stable numerically than the results of the standard 
analysis in the MS-scheme. Then numerical values for 
the reference MS-scheme parameters can be obtained by 
a renormalization group "rotation" from the numerical 
values found in the internal schemes. Renormalization 
group "rotation" (the re-calculation of numerical values 
from one scheme to another) is a quite formal operation 
and can be easily controlled numerically. One example 
of such a "rotation" (the renormalization group scaling 
which is a one-parameter subgroup of the general renor- 
malization group) is the evolution of the coupling con- 
stant to the reference scale Mz- Below we give a detailed 
description of our approach. 

The normalized r-lepton decay rate into nonstrange 
hadrons hs=o is given by 



cq. (||) is the parton model result while the second term 
Sp represents perturbative QCD effects. For the flavor 
mixing matrix element we use |K»z| 2 = 0.9511 ± 0.0014 
||. The factor Sew — 1-0194 is an electroweak cor- 
rection term |28| and Sew = 0.001 is an additive elec- 
troweak correction P^| . The nonperturbative correc- 
tions are rather small and consistent with zero; we use 
S NP = -0.003 ±0.003 (see e.g. g|). Note that recently 
the problem of duality violation for two-point correla- 
tors has been discussed However, no established 
quantitative estimates of that violation are available yet. 
Considerations show that they can be rather large and 
can reach the level of few percents. This problem can 
affect the numerical value of the coupling extracted from 
the analysis because of the numerical change of the quan- 
tity Sp extracted from eq. (||). In the present note we 
concentrate on the perturbative part of the decay rate 
and numerical uncertainties related to the renormaliza- 
tion scheme freedom of perturbation theory. In this re- 
spect new possible corrections do not qualitatively affect 
our analysis. The corrections due to duality violation are 
of a new nature and they can be added independently to 
eq. (|j|). They would only change the input numerical 
value for the Sp within our approach. 

The value for the decay rate R T s=o has been measured 
by the ALEPH and OPAL j| collaborations with re- 
sults very close to each other. For definiteness we use 
the ALEPH data and briefly comment on the OPAL data 
later on. With the experimental result 



R e T x s p =0 = 3.492 ± 0.016 



one obtains from eq. (^) 



Sp p = 0.203 ± 0.007 . 



(4) 



(5) 



The basic object of the theoretical calculation is Adler's 
-D-function which is computable in perturbation theory 
in the Euclidean domain. In the MS-scheme the pertur- 
bative expansion for the D-function is given by 



+k 3 



a s (Q) 



0(a s (Q) 5 ) 



(6) 



with (see e.g. U) 
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*i = ^-9C(3), 

58057 779 75 

fc2 = W-- c(3) + T c(5) - 



(7) 



_ r(r -» h s= ov) 
>n " T(t->Iuv) 

= N c \V ud \ 2 S E w(l + Sp + S E w + S N p) (3) 



where N r = 3 is the number of colors. The first term in 



Here £(x) is Riemann's ^-function. I n the following we 
use the notation 



a„{Q) 



a s {Q) 



(8) 
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for the standard MS-coupling constant normalized at the 
scale jj. = Q. Numerically we find 



D{Q 2 ) = 1 + a s {Q) + 1.6398a s (Q) 2 
+k 3 a s (Q) 4 + 0{a 5 s {Q)) . 



6.3710a s (Q) E 



(9) 



The coefficient k% is still unknown which prevents us from 
using the last term in eq. (^J) for our analysis. We nev- 
ertheless list this term throughout the paper to obtain a 
feeling for the possible magnitude of the 0{a A s ) correc- 
tion. The particular numerical value of k% ~ 25 is ob- 
tained on the basis of geometric scries approximation for 
the series @ and is often used in the literature |l6|,|32|,^3| . 
In our analysis we do not use any particular numerical 
value for k 3 and only give some illustrative results of the 
influence of this term on the numerical value of the cou- 
pling constant extracted from r-data. 

In the MS-scheme the perturbative correction Sp is 
given by the perturbation theory expansion 



5.2023< + 26.366< 



+(78.003 + fc 3 K + 0(a 5 s 



(10) 



where the MS-scheme coupling constant a s = 7ra s is 
taken at the scale of the T-lepton mass fj, = M T — 
1.777 GeV. Usually one extracts a numerical value for 
a s (M T ) by treating the first three terms of the expression 
in eq. ( |lO|) as an exact function - the cubic polynomial, 
i.e. one solves the equation 



a s + 5.2023a 2 
The solution reads 



26.366a 3 ; = S e P xp . 



iraf (M T ) = af{M T ) = 0.3404 ± 0.0073 



exp • 



(ii) 



(12) 



We call this method the standard method. The quoted 
error is due to the error in the input value of 5 e p xp . We 
retain some additional decimal points in the numerical 
expression for the coupling constant in order to use them 
for the evolution of the coupling constant to the scale 
Mz- It is rather difficult to estimate the theoretical un- 
certainty of the procedure itself. The main problem is to 
estimate the quality of the approximation for the (asymp- 
totic) series in eq. ( |To| ) given by the cubic polynomial in 
eq. (0). 

As a criterion of the quality of the approximation one 
can use the pattern of convergence of the series (|l0]) which 
is 



S e p xp = 0.203 = 0.108 



0.061 + 0.034 



(13) 



One sees that the corrections provide a 100% change of 
the leading term. Another criterion is the order-by-order 
behavior of the extracted numerical value for the coupling 
constant. In consecutive orders of perturbation theory 
(LO - leading order, NLO - next-to-leading order, NNLO 
- next-next-to-leading order) one has 



af(M T ) LO = 0.6377, af{M T ) NLO 
a s *{M T ) NNLO =0.3404. 



0.3882, 



(14) 



Formally we obtain a series for the numerical value of the 
coupling constant of the form 



af(M T ; 



NNLO 



= 0.6377- 0.2495 - 0.0478 



(15) 



Limiting ourselves to the next-to-next-to-leading order 
result (NNLO) we can take a half of the last term as an 
estimate of the theoretical uncertainty. It is only an in- 
dicative estimate. No rigorous justification can be given 
for such an assumption about the accuracy of the ap- 
proximation without knowledge of the structure of the 
whole series. Nevertheless we stick to this definition for 
our purposes. The theoretical uncertainty obtained in 
such a way - Aaf (M T ) ( ft = 0.0239 - is much larger than 
the experimental uncertainty given in eq. (12). This is 
a challenge for the theory: the accuracy of theoretical 
formulas cannot compete with experimental precision at 
present. Assuming this theoretical uncertainty we have 



a** (Mt)nnlo — 

0.3404 ± 0.0239 t/l ± 0.0073 e 



(16) 



Theory dominates the error even if the estimate for its 
precision ±0.0239^ is not reliable (heuristic and only in- 
dicative). Thus the straightforward analysis in the MS- 
scheme is not stable numerically and the naive estimate 
of the theoretical uncertainty is large. 

The use of the MS-scheme is not obligatory for practi- 
cal calculations. The MS-scheme has a history of success 
for massless calculations where its results look natural 
and the corrections are usually small. This is not the 
strict rule, however, and there are cases (like gluonic cor- 
relators |54|) where corrections dramatically depend on 
the quantum numbers of the operators. In fact, the MS- 
scheme is rather artificial. It is simply defined by con- 
vention (let us be remindful of the evolution from the 
MS-scheme to the MS-scheme which had its origin only 
in technical convenience [[35]]). From technical point of 
view, in practical calculations of massless diagrams of 
the propagator type, another scheme - the G-scheme - is 
the most natural one Q . It normalizes the basic quan- 
tity of the whole calculation within integration-by-parts 
technique - one loop masless scalar diagram - to unity 
p7| . /3-functions coincide in both schemes. It could have 
well happened that the G-scheme would be historically 
adopted as the reference scheme because corrections in 
this scheme are typically smaller than that in the MS- 
scheme. However, for the tau system the direct (stan- 
dard) analysis in the G-scheme fails. 

Therefore different schemes used for the numerical 
analysis can produce rather different numerical results 
for the final reference quantity - the coupling constant in 
the MS-scheme. Note that strictly speaking any scheme 
is suitable for a given perturbative calculation. However, 
it can lead to unusual (or even unacceptable) results in 
a numerical analysis. The only criterion for the choice 
of scheme at present is the heuristic requirement of fast 
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explicitl convergence: the terms of the series should de- 
crease. Clearly this is a rather unreliable criterion. It 
does not provide strict quantitative constraints necessary 
for the level of precision usually claimed for the r-system 
analysis. 

In the following we suggest a new procedure for ex- 
tracting a s in the MS-scheme from the r system without 
explicit use of eq. (10). This procedure is applicable to 
any observable in whatever scheme it was originally com- 
puted. The observation is that any perturbation theory 
observable generates a scale due to dimensional trans- 
mutation and this is its internal scale. It is natural for a 
numerical analysis (and is our suggestion) to determine 
this scale fisrt and then to transform the result into a MS- 
scheme parameter (or any other reference scheme) using 
the renormalization group invariance. We deliberately 
use the explicit renormalization scheme invariance of the 
theory to bring the result of the perturbation theory cal- 
culation into a special scheme first, then we perform a 
numerical analysis in this particular scheme. Only after 
that we transform the obtained numbers into the refer- 
ence MS-scheme. The last step is done only for compar- 
ison with other experiments (or just for convenience; the 
system itself can be well described in its internal scheme 
without any reference to the MS-scheme). This is our 
suggestion for the resolution of the problem of numeri- 
cal instability of extracting parameters from truncated 
perturbation theory expressions. 

A dimensional scale in QCD emerges as a boundary 
value parameterizing the evolution trajectoty of the cou- 
pling constant. The renormalization group equation 



^ 2 3~2 a (^ 2 ) = P{ a {v 2 )) , a=- 



is solved by the integral 

'"(f)=*«" 2 » + / 



/ 1 



1 



(17) 



di (18) 



where the indefinite integral $(a) is normalized as follows 



*(a) 



1 , r 1 ,fh } 
/MO aA) $ n VA> + a/?i 



a(3 2 



(19) 



Here /^(a) and j3(a) denote the second order and full (3 
function, or as many terms as are available, given by 

/3 2 (a) = -a 2 (p +aPi), 

(3(a) = -a 2 (f3o + f3 ia + (3 2 a 2 + f3 3 a 3 ) + 0(a 6 ) , (20) 

a is a generic coupling constant. The four- loop /3-function 
coefficient j3 3 is now known in the MS-scheme [B8|| 



Ih 



140599 445 

1 C 3) = 47.228. . 

4608 32 ^ v ; 



(21) 



The integration constant in eq. (18) is adjusted such that 
the asymptotic expansion of the coupling constant at 
large momenta Q 2 — > oo reads 



a(Q 2 ) = 



L = In 



PoL [ l (3 2 
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A 2 



L 2 



O 



L 3 



(22) 



This serves to define the parameter A (dimensional scale) 
for a generic coupling constant. A s is the standard MS- 
scheme scale for the coupling constant a s . 

The solution ( |T8| ) of the renormalization group equa- 
tion (|l7| ) describes the evolution trajectory of the cou- 
pling constant. This trajectory is parametrized by the 
scale parameter A and the coefficients of the (3 function 
j3i with i > 2 (see e.g. |39|). The evolution is invariant 
under the renormalization group transformation 

a — > a(l + Ki<j + n 2 a 2 + n 3 a 3 + . . .) (23) 

with the simultaneous change 

A 2 -> A 2 e~ Kl / /3 ° , (24) 

/?o,i left invariant and 

I3 2 — ► Aj - kI/3 + k 2 /?o - KlA 

-»■ #3 + 4k?/3 + 2k 3 /3 + k\(3 x - 2 Ki (3k 2 /3 +&). 

If this transformation was considered to be exact and 
the exact /3-function corresponding to the new charge 
was used then it would be just a change of variable in a 
differential equation dl7|) or the exact reparametrization 
of the trajectory ( |l8| ) and hence would lead to identical 
results. However, the renormalization group invariance 
of eq. ([l^) is violated in higher orders of the coupling 
constant because we consistently omit higher orders in 
the perturbation theory expressions for the /3-functions. 
This is the point where the finite-order perturbation the- 
ory approximation for the respective /^-functions is made. 
This is the source for different numerical outputs of anal- 
yses in different schemes. 

Our procedure for the extraction of a s is heavily based 
on the formal renormalization group invariance of the 
theory. We claim that because of this invariance we can 
do our numerical analysis in any scheme. The reason 
for the choice of a particular scheme is only the quality 
of the convergence (which, of course, is subject to some 
personal taste). We have chosen the effective scheme 
because we consider it to be more consistent and more 
stable numerically. 

Technically we introduce an effective charge a T through 
the relation (27|4(]-{4§ 



5 f =a T = ^ 

7T 



(25) 



and extract the parameter A T which is associated with a T 
through eq. (18). This is just the internal scale associated 
with the physical observable R T . The effective (3- function 
is given by the expression 

(3 T = -a 2 T {/3 T0 + (3 Tl a T + (3 t2 cl 2 t + /3 r3 a 3 T + . . .) (26) 
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with p r0 = fa, Pti = Pi, and 



af SI (M T ) = 0.3184 ±0.0060 e 



(29) 



Pt2 



-12.3204, p r3 



•182.719+ -k 3 . (27) 



The extraction of the numerical value for the internal 
scale A r is done from equation ( |rs| ) with a T (M T ) = 8 £ p P . 
The coefficient p r3 does not enter the analysis. The pa- 
rameter A s = Aj^g is found according to eq. (p4|). The 
MS coupling at fi = M T is obtained by solving eq. (O) for 
a s (M T ) with regard to ln(M^/A|) which is known if A s 
is obtained; the /3-function is taken in the MS-scheme. 
For consistency reasons we only use the MS-scheme /?- 
function to three-loop order since the effective /3-function 
P r is only known up to the second order, cf. eq. (27). 
A N 3 LO analysis is possible only if a definite value is 
chosen for k 3 . We give some estimates later. 

Our procedure is based on renormalization group in- 
variance and one can start from the expression for the 
decay rate obtained in any scheme. The only pertur- 
bative objects present are the /3-functions. Both p^ 
and P T , however, converge reasonably well which is the 
only perturbation theory restriction in our method. It 
also highlights the limit of precision within our proce- 
dure: the expansion for p T is believed to be asymptotic 
as any expansion in perturbation theory. The asymp- 
totic expansion provides only limited accuracy for any 
given numerical value of the expansion parameter which 
cannot be further improved by including higher order 
terms. The expansion used is presumably rather close to 
its asymptotic limit as can be seen by taking a look at 
the expansion 



/9 

p T (a T ) = -<4(-7 +4a r - 12.3204(4 



! -182.719+ ~k 3 ) ) + 0(<) 



(28) 



with a T ~ 0.2 at the scale M T . The convergence of the 
series depends crucially on the numerical value of £3. If 
k 3 had a value where the asymptotic growth starts at 
third order then further improvement of the accuracy 
within finite-order perturbation theory is impossible. 

At every order of the analysis we use the whole in- 
formation of the perturbation theory calculation. Espe- 
cially, the appropriate coefficient of the /3 T -function is 
present. In the standard method the coefficient P2 enters 
only at order 0(a 4 s ) of the r-lepton decay rate expansion. 
We call our procedure the renormalization scheme invari- 
ant extraction method (RSI) hoping that it is clear what 
is meant by this name from our explanations. Note also 
that a s itself is not a physical object and is renormaliza- 
tion scheme noninvariant. In this respect we extract the 
noninvariant parameter a s using invariance of the physics 
in order to perform the numerical analysis in the most 
suitable scheme. Then the output of the analysis is sim- 
ply transformed into a numerical value for a s according 
to the renormalization group transformation rules. For 
the coupling constant in the MS-scheme in NNLO we find 



which is smaller than the corresponding value obtained 
within the standard procedure eq. ([[2]) . How to estimate 
the quality of this result? The parameter which is really 
extracted in consecutive orders of perturbation theory 
within our method is the scale A T . Because of the relation 
(see eqs. 

A s = A T e" 5 - 20232/2/3 ° = 0.3147A T (30) 

we can look at A s directly. We find 

A s \lo = 595 MeV, A s \nlo = 288 MeV, 
A s \ NN lo = 349 Me V (31) 

or, representing the NNLO result as a formal series, 

A s \nnlo = 595 - 307 + 61 - . . . MeV. (32) 

Note that at leading order the scales (as well as charges) 
are equal in all schemes. Therefore the leading order 
result (A s \lo = 595 MeV) is not representative, only 
indicative. Assuming according to our convention that 
the uncertainty of A s is given by the half of the last term 
of the series ( p2| ) we have 



A, = 349 ± 31 MeV 



(33) 



which leads to the numerical value for the MS-scheme 
coupling constant 



°- 3184 +aoi60 ■ 



(34) 



This result is obtained from eq. (|18|) with three-loop /3- 
function. Taking the average we find 



0.3184 + 0.0159. 



(35) 



This is better than the theoretical error of the standard 
result eq. (|l6|). Still the theoretical error should be con- 
sidered as a guess rather than a well-justified estimate of 
the uncertainty. 

Let us briefly comment on the k 3 contribution. Clearly 
the estimate k 3 = 25 is rather speculative. We, therefore, 
use a different strategy in the analysis. We determine 
the range of k 3 which is safe for explicit convergence of 
perturbation theory. If the actual value of k 3 will be 
discovered in this range then perturbation theory is still 
valid and will give better accuracy in NNNLO. If not, 
the asymptotic growth of perturbation theory series is 
already reached and its accuracy cannot be improved. 

We require that the last term is equal to the half of 
the previous one. In the standard way (eq. [HI) we have 



|(78 + fc 3 )a s | < - 26.36 « 13 



which for a s =0.1 gives 



- 208 < fcf* < 52 . 



(36) 



(37) 
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In the RSI way (eq. |2^) we have 

|(-182+^fc 3 )a T | < 1 12.32 « 6 



which for a T = 0.2 gives 



33. 



<kl< 47.1 



This range is much narrower than that in eq. (|37|). The 
effective scheme method is much more sensitive to the 
structure of the series as can be seen from eq. (p8|). The 
actual precision depends on the actual value chosen for ks 
and it is rather premature to speculate about numbers. 

Still we show the worst result (in the optimistic sce- 
nario that &3 lies in the safe range) that can be expected 
within the RSI approach. In the RSI approach with 
k% = 47 we find the scale parameter in NNNLO 



rif =4 effective theories is done directly at the scale M T , 
which is possible, then the result is slighly smaller than 

(38) in the case when the evolution within nf = 3 effective 
theory is done first to the scale \i c . In the following we 
stick to the procedure where the matching is performed 
precisely at the matching scales fj, c , b - We first run the 

(39) coupling constant within nf = 3 effective theory from 
the scale M T to fi c then match the result to nf = 4 
coupling constant, run it to \x b and match to nf = 5 cou- 
pling constant. The last step is just evolution to Mz- 
Note that the alternative would be to perform matching 
between ny = 3 and nj — 4 effective theories directly at 
the scale M T (because it is rather close to a c ) but in this 
case the final result is slightly smaller than in our present 
procedure. 

The running to the scale Mz gives the following result 
for the standard method estimate 



A 



s I N N N LO 



334 MeV . 



(40) af(M z ) = 0.1210 ± 0.0008 ea;p ± 0.0006 c ± 0.0001 6 (45) 



With &3 = 34 one has 

A s \nnnlo = 367 MeV. 
Taking the average we have 

A, = 350 ± 17 MeV 



(41) 



(42) 



which is the best possible estimate if we require that the 
perturbation theory series for the /3 r -function still con- 
verges (according to our quantitative criterion of con- 
vergence). That results in the numerical value for the 
MS-scheme coupling constant found with four-loop /3- 
function from eq. (|18| ) 



0.3133 < a, < 0.3314. 



(43) 



Therefore our conservative estimate of the theoretical er- 
ror in the optimistic scenario for the convergence of per- 
turbation theory series in NNNLO reads 



a, = 0.322 ±0.009. 



(44) 



While the estimation of the theoretical uncertainty is a 
tricky matter and can be considered as indicative the cen- 
tral numerical value of the coupling constant definitely 
becomes smaller as compared to the standard result. 

At present the reference value for the coupling con- 
stant is commonly given at the scale Mz — 91.187 GeV. 
The running to this reference scale is done with the 



four-loop f3- function in the MS-scheme |38 and three- 
loop matching conditions at the heavy quark (charm and 
bottom) thresholds [Q. For the threshold parameters 
related to heavy quark masses we use u c = fh c (fj, c ) — 
(1.35 ± 0.15) GeV and u b = rh b (u b ) = (4.21 ± 0.11) GeV 
(e.g. jf5|) where rh q (fi) is the running mass of the heavy 
quark in the MS-scheme. Note that because of the trun- 
cation of matching conditions the result of the running 
slightly depends on at what scale the matching is actu- 
ally performed. If the matching between the nf = 3 and 



exp 

where the subscript exp denotes the error originating 
from 8 e p V . The errors with subscripts c, b arise from 
the uncertainty of the numerical values of the charm and 
bottom quark masses that enter the evolution analysis. 
These errors are rather small (we retain the additional 
decimal place in the result, which is not really justified 
from the precision of the experimental input, just to show 
these uncertainties). If the matching between the nf = 3 
and nf = 4 effective theories is done directly at the scale 
M T one has to change the central value 0.1210 — > 0.1202 
which shows the uncertainty related to the truncation of 
the matching conditions. 

The central value in eq. ( fi"5| ) is slightly higher than that 
calculated from high energy experiments [|3|. The theo- 
retical perturbative expansions for observables in high 
energy experiments converge better numerically than ex- 
pansions at low energies because the coupling, which is 
the parameter of the perturbative expansion, is smaller 
at higher energies due to the property of asymptotic free- 
dom. This feature makes it less important to treat the 
higher order terms carefully in high energy applications 
as compared to the low energy r-lepton estimates. How- 
ever, the experimental data in high energy experiments 
are usually less precise which leads to large errors in the 
a s determination from high energy experiments. The 
fact that the value in eq. (|4^) is higher than that calcu- 
lated from high energy experiments caused some discus- 
sion about the reliability of estimates from the r-lepton 
data. Our analysis resolves this problem. The running 



of (M T ) given in eq. ( |35| ) to Mz with the four-loop 
/3-function and with three-loop heavy quark matching ac- 
curacy gives 



a 



RSI 



(M z ) = 0.1184 ±0.00074 



exp 



±0.00053 c ± 0.00005& 



(46) 



where we have kept five decimal places in order to ex- 
hibit the magnitude of different sources of uncertainty. 
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Eq. ( }46| ) constitutes our main result for the coupling 
a s {Mz) derived from tau data. 

The OPAL collaboration has reported an experimental 



value of R e T x s p =0 = 3.484 ± 0.024 ||. This leads to 8p x p = 
0.200 ± 0.009 eKp and 



af SI (M T ) = 0.3158 ± 0.0078 ea: p 

which, when evolved to Mz, gives 

af SI (M z ) = 0.1181 ± 0.00097 ea; p 

±0.00052 c ± 0.00005 6 . 



(47) 



(48) 



This value is close to the one in eq. (|46|) based on the 
ALEPH data. 

The theoretical uncertainty comes mainly from the 
truncation of the perturbation theory series. Taking the 
result of the NNLO analysis eq. (pq) we find 



A, 



..RSI 



(Mz)th = 0.0019 



(49) 



In the most optimistic scenario with the NNNLO analysis 
eq. (Eh one has 



af SI (M z ) N 3 LO = 0.119 ± 0.001 



(50) 



As we have already noted the interpretation of the 
higher order terms in the perturbation theory expansion 
is numerically important for the analysis of the r-data. 
The regular method to resum higher order perturbation 
theory corrections is based on the direct integration of 
the renormalization group improved correlators over the 
contour in the complex Q 2 plane |l4| . This method al- 
lows one to resum corrections generated by the running 
of the coupling constant along the integration contour 
and is now widely used for the analysis of the r-data. 
We now briefly comment on the extraction of the strong 
coupling constant within rcsummed perturbation theory. 
As in ref. |^6| we fit the theoretical expression for the 
decay rate in the contour improved approach to the ex- 
perimental result 8p xp eq. (|^) and find 



a.y{M T ) = 0.343 ±0.009 e 



(51) 



within the renormalization scheme invariant extraction 
method described above i.e. with the introduction of the 
effective charge first. This value differs from the finite- 
order perturbation theory result eq. (p9|). Note that the 
two values extracted from finite-order perturbation the- 
ory analysis eq. ( p9| ) and the contour improved perturba- 
tion theory analysis eq. (51) do not overlap within their 
respective error bars given from the experimental uncer- 
tainty only. This situation was anticipated in [ fl4| where 
the resummed NNLO analysis had been first performed. 
The point is clear: resummation provides a specific esti- 
mate of higher order terms. In finite-order perturbation 
theory one adopts a model where all higher order terms 
have been neglected. In contour improved perturbation 
theory one adopts an explicit model with higher order 



terms generated by the running of the coupling constant 
along the integration contour. With present experimen- 
tal accuracy one can already distinguish between these 
two possibilities. One should always keep in mind that 
the two determinations eq. ( p9| ) and eq. ( |5l| ) result from 
different models and one should not mix their predictions. 

The numerical value of the coupling constant appropri- 
ate for high energy experiments is normally small (much 
smaller than for r-data) and perturbation theory con- 
verges faster (in similar kinematical situations). The re- 
summation does not produce any big numerical changes. 
Therefore finite-order perturbation theory is normally 
used for the analysis of high energy experiments (resum- 
mation of the contour type can be done but produces a 
small numerical effect) and one usually quotes numeri- 
cal values of the coupling constant extracted with finite- 
order perturbation theory. Or resummation of the sort 
different from that used for the r system is used (like 
Coulomb type resummation for heavy quarks f47j , |48|| ) . 
Therefore we suggest to use the finite-order perturba- 
tion theory prediction for the coupling constant extracted 
from r-data in order to compare it with the results of high 
energy experiments. 

To conclude, we have extracted the numerical value of 
the strong coupling constant from r-data within a proce- 
dure based on explicit use of renormalization scheme in- 
variance. The numerical value for the coupling constant 
is systematically smaller than that derived by the stan- 
dard treatment. When evolved to Mz our MS-scheme 
value for the coupling constant extracted in finite-order 
perturbation theory reads 

a s (M z ) = 0.1184 ± 0.0007 ea; p ± 0.0006/* mass ■ (52) 

This central value is closer to the value of a s derived from 
high energy experiments than previous determinations of 
a s from r-data. The theoretical uncertainty of the result 
is still only indicative: it ranges from the conservative es- 
timate in NNLO Aa s (Mz)th = ±0.0019 to an optimistic 
one based on the assumption about NNNLO contribution 
Aa s (M z )th = ±0.001. 
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